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INNER  PRODUCT  PERFORMANCE  CRITERIA  FOR  EVALUATING 
COMBAT  MODELS 

HERBERT  E.  COHEN 


U.S.  ARMY  MATERIEL  SYSTEMS  ANALYSIS  ACTIVITY 
ABERDEEN  PROVING  GROUND,  MARYLAND  21005-5071 


The  U.S.  Army  is  actively  engaged  in  developing  combat  models  using 
approximate  solutions  to  the  linear  homogeneous  Lanchester  equations.  How 
well  these  approximate  solutions  deviate  from  the  true  or  exact  solution  is  an 
important  issue  within  the  Arm>  if  we  are  going  to  have  any  confidence  in 
these  approximate  solutions.  An  inner  product  performance  criteria  is 
proposed  which  has  a  simple  geometric  meaning  and  is  easy  to  calculate  for 
evaluating  how  well  the  approximate  solution  agrees  with  the  exact  solution 
associated  with  the  linear  homogeneous  Lanchester  equation  of  combat. 

It  is  shown  that  an  inner  product  performance  criteria  is  a  useful  measure  to 
evaluate  the  performance  of  current  combat  models  against  the  exact  solutions 
for  piecewise  linear  homogeneous  Lanchester  equations  for  all  homogeneous  or 
heterogeneous  battles. 

The  Lanchester  homogeneous  equations  in  vector  form  may  be  given  by  Eq  (1) 

,v(f)-^(0  *(*)  (1) 

where  x  is  an  n  x  1  vector  such  that 

,y(0  =  (  R  !  R  2. .  .R  p,  B  !  ,  B  2> 

Ri  is  the  i^  Red  weapon  system  strength  at  time  t. 

Bj  is  the  jth  Blue  weapon  system  strength  at  time  t. 
and 

p  +  q  =  n 

Ri  and  Bj  both  being  nonnegative, 

with  an  initial  condition  x(t0)  =  x0,  A(t)  is  the  time  varying  n  by  n 
attrition  coefficient  matrix.  Assuming  A(t)  is  piecpwise  constant  for 
t i <t< t i + i  for  i  =0,  1  ...,  n  where  tn+j  =  tf,  then  Eq  (1)  has  an  exact 
solution  given  by 


•v(/,*,)  =  f>  A  •(/,)  (2) 

for  i  =  0 ,  l  ,  . . . ,  a? 

We  now  define  the  adjoint  differential  equation  of  Equation  (1)  given  by 

y  =  -i'(/)y  (3) 

with  a  value  of  y(t)  at  t  =  tf  given  by  y  (tf). 

We  now  claim  that  the  inner  product  yT(t)  x  (t)  is  a  constant.  This  can  be 
shown  as  follows: 

^yT(0.v(0  =  yT-v(0  +  yT(0.v(0 

=  ( -  .d 7  y ) r  -v  ( o  +  y 7  ( O  a*  ( 0 

=  y7(-/Lv)  +  yT(0.v(0 
=  -  y 7  ( O  v  ( o  -  y 7  ( O  v  ( o 
=  0  QZTD. 

where  we  have  used  Eq  (1)  and  Eq  (3). 

Thus 

y  (/)  v ( / )  =  constant  (4) 

or  more  specifically  at  the  final  time,  tf,  of  the  battle 

yT(//).v(//)  =  yr(0)x(0)  (5) 
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In  fact  if  we  select  y(tf)  to  be  orthogonal  to  x  (tf)  we  have 


7 '  (//)■'■(//)  =  0  (6) 

and  the  inner  product  yT(t)  x  (t)  will  remain  zero  for  all  t.  Hence  all 
solutions  of  the  adjoint  equation,  y(t)  in  Eq  (3),  are  orthogonal  to  the 
surface  containing  the  solution  of  the  Lanchester  equation  in  Eq  (1),  even  for 
time  varying  attrition  matrix  and  all  battles.  The  geometric  meaning  is 
simply  that  the  projection  of  x(t)  on  the  adjoint  solution  y(t)  multiplied  by 
the  magnitude  of  y(t)  is  zero.  Figure  1  provides  the  geometric  relationship 
and  the  position  of  a  generic  Army  Model  solution,  x/\^,  in  relationship  to  x 
and  y . 

Within  the  Army,  the  model  given  by  x/\^,  for  Army  Model,  uses  the  Lanchester 
equation  to  describe  many  combat  situations.  The  finite  difference  equation 
for  Eq  (1)  is 

■V  ,u  (/  +  . AO -(/  +  . -I -\/).v 


u  hich  u  ilh  A  /  =  1 


V  |  +  “  (/  +  1  )  X  AM  (Y) 


If  is  to  be  a  close  approximation  to  x,  then  the  inner  product  of  x/\\j\  with 
y  of  the  adjoint  equation,  Eq  (3),  is  the  projection  of  x/\m  along  y  and  thus 
the  distance  from  the  true  solution. 

Data  is  presented  in  Table  1,  x/\m  is  in  the  second  column  and  the  exact 
solution,  x(t),  generated  by  Eq  (2)  over  each  time  interval  is  presented  in 
Figures  2  to  6.  Runge  Kutta  solutions  are  also  presented  in  these  figures. 

Table  2  provides  the  results  of  the  inner  product  y^(t)x(t)  and  yT(t)  x/\m  over 
each  time  interval.  Clearly  y^(t)  x/\m  is  not  constant  and  progressively  moves 
away  from  the  exact  solution  x(t)  with  time.  This  is  shown  in  Figure  7.  We 
have  a  clear  measure  of  how  x/\m  deviates  from  the  true  solution.  Thus,  the 
inner  product  criteria  offers  a  simple  method  for  testing  any  linear 
homogeneous  Lanchester  model  of  combat  for  all  battles. 

In  conclusion,  we  have  demonstrated  that  the  inner  product  criteria  using  the 
adjoint  differential  equation  can  provide  combat  modelers  with  a  useful  simple 
criteria  for  evaluating  the  performance  of  approximate  solutions  to  the  linear 
homogeneous  Lanchester  equation  of  combat. 
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FIGURE  1 


X  =  EXACT  LANCHESTER  SOLUTION 
Y  =  ADJOINT  SOLUTION 
XAM=  ARMY  MODEL  SOLUTION 

YtXam=  MEASURE  OF  DISTANCE  FROM  TRUE  SOLUTION  X. 
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TABLE  1.  DATA. 


Two-on-Three  Scenario 

(With  Suppression  Off  and  100%  in  Quadrant) 

1  Blue  Initial  Strength,  50 

2  Blue  Initial  Strength,  50 

3  Red  Initial  Strength,  50 

4  Red  Initial  Strength,  25 

5  Red  Initial  Strength,  25 


No.  of 

Tarqet 

F  i  rer 

F i rers 

1 

2 

3 

4 

5 

Time  =  0  minutes 

IB 

50.0000 

0. 

0. 

0.0063643 

0.0000000 

0 . 0  0  0 1 U 1 4 

2B 

50.0000 

0. 

0. 

0.0024091 

0. 

0. 

3R 

50.0000 

0.0083374 

0. 

0. 

0. 

0. 

4R 

25.0000 

0.0008168 

0. 

0. 

0. 

0. 

5R 

25.0000 

0. 

0. 

0. 

0. 

0. 

Time  =  1  minute 

IB 

23.7626 

0. 

0. 

0.0043994 

0.0000000 

0.0023342 

2B 

50.0000 

0. 

0. 

0.0023752 

0.0000000 

0.0000519 

3R 

23.6800 

0.0083374 

0. 

0. 

0. 

0. 

4R 

0 . 0C0D1C8 

0. 

n. 

0. 

0. 

5R 

24.6958 

0. 

0. 

6. 

0. 

0. 

Time  =  2  minutes 

IB 

10.6916 

0. 

0. 

0.0024300 

0.0000816 

0.0048297 

2B 

50.0000 

0. 

0. 

0.0020609 

0.0000124 

0.0005698 

3R 

10.2977 

0.0083343 

0.0000027 

0. 

0. 

0. 

4R 

24.9896 

0.0008168 

0. 

0. 

0. 

0. 

5R 

21.2119 

0. 

0. 

0. 

0. 

0. 

Time  =  3  minutes 

IB 

4.3175 

0. 

0. 

0.0004420 

0.0005786 

0.0068515 

2B 

49.9983 

0. 

0. 

0.0007079 

0.0001323 

0.0024820 

3R 

2.6462 

0.0079945 

0.0003013 

0. 

0. 

0. 

4R 

24.9000 

0.0008046 

0.0000623 

0. 

0. 

0. 

5R 

16.4044 

0. 

0. 

0. 

0. 

0. 

Time  =  4  minutes 

IB 

1.8461 

0. 

0. 

0.0000143 

0.0035953 

0.0034569 

2B 

49.8574 

0. 

0. 

0.0001091 

0.0009374 

0.0025676 

3R 

0.4079 

0.0062070 

0.0018716 

0. 

0. 

0. 

4R 

24.3531 

0.0006316 

0.0006921 

0. 

0. 

0. 

5R 

7.1839 

0. 

0. 

0. 

0. 

0. 

Red  unit  removed  after  4  minutes. 
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Firer  1  obi 
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Firer  3  ori 


o 

o 

o 

o 

o 

O 

o 

CO 

10 

CO 

CM 

HioNaais 

3  DHOd 

8 

Series  1  ~+~  Series  2  Series  3 

RUNGE-KUTTA  EXACT  SOLUTION  ARMY  MODEL 


Firer  5  or) 
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Series  1  Series  2  Series  3 

runge-kutta  exact  solution  army  model 


TABLE  2. 


yT(t)x(t) 

t 

Initial  t=0 

Initial  t=5 

0 

0 

-4.16  x  10'15 

1 

-2.58  x  10'14 

1.48  x  10'15 

2 

-1.914  x  10‘14 

-4.07  x  10*15 

3 

-3.689  x  10'14 

-6.31  x  10'17 

4 

-4.893  x  10'14 

2.49  x  10'16 

yT(t)  XAM 

t  yT 

(t)  XAM 

0 

0 

1 

-8.6 

2 

-14.26 

3 

-17.756 

4 

1  -18.7 
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SUBJECT:  Inner  Product  Performance  Criteria  for  Evaluationing 
Combat  Models. 

REASON  FOR  PERFORMING  THIS  EFFORT:  Investigation  of  the  properties  of  the 
linear  homogeneous  Lanchester  equations  which  are  used  extensively  by  Army 
combat  modelers. 

MAIN  OBJECTIVES  OF  THE  EFFORT:  Investigate  invariant  properties  of  liniar 
homogeneous  Lanchester  equations  of  combat  and  to  compare  theoretical  results 
with  those  derived  from  Army  combat  models. 

SCOPE  OF  EFFORT:  Develop  solutions  for  piecewise  constant  linear  homogeneous 
Lanchester  equations  in  state  variables  and  introduce  the  role  of  the 
adjoint  differential  equation. 

IMPACT  OF  THE  EFFORT:  Provides  explicit  solution,  independent  of  step  size, 
of  the  linear  homogeneous  lanchester  equations  thus  increasing  speed  of 
computation  and  develops  a  simple  performance  criteria  for  evaluating  how 
well  the  approximate  solutions  agree  with  the  exact  solutions  of  the 
equation  of  combat. 
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PERFORMING  ORGANIZATION/POINT  OF  CONTACT:  U.S.  Army  Materiel  Systems 
Analysis  Activity  (AMSAA) ,  Special  Studies  and  Activities  Office, 

Herbert  E.  Cohen,  DSN  298-2785. 
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